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Abstract 

Golay sequences are well suited for use as codewords in orthogonal frequency-division multiplexing 
(OFDM) since their peak-to-mean envelope power ratio (PMEPR) in g-ary phase-shift keying (PSK) 
modulation is at most 2. It is known that a family of polyphase Golay sequences of length 2™ 
organizes in m!/2 cosets of a q-ary generalization of the first-order Reed-Muller code, RMg(l,m). 
In this paper a more general construction technique for cosets of RMq(l,r7i) with low PMEPR is 
established. These cosets contain so-called near-complementary sequences. The application of this theory 
is then illustrated by providing some construction examples. First, it is shown that the m!/2 cosets of 
RMg(l, m) comprised of Golay sequences just arise as a special case. Second, further families of cosets 
of RMq(l,m) with maximum PMEPR between 2 and 4 are presented, showing that some previously 
unexplained phenomena can now be understood within a unified framework. A lower bound on the 
PMEPR of cosets of RMq(l, m) is proved as well, and it is demonstrated that the upper bound on the 
PMEPR is tight in many cases. Finally it is shown that all upper bounds on the PMEPR of cosets of 
RMg(l, m) also hold for the peak-to-average power ratio (PAPR) under the Walsh-Hadamard transform. 

Index Terms 

Aperiodic, code, complementary, correlation, Golay, OFDM, PAPR, PMEPR, Reed-Muller, Rudin- 
Shapiro, sequence, Walsh-Hadamard 

I. Introduction 

Despite many evident advantages of the orthogonal frequency-division multiplexing (OFDM) 
modulation technique, the widespread acceptance of OFDM mainly suffers from the usually 
high peak- to-mean envelope power ratio (PMEPR) of uncoded OFDM signals. A promising and 
elegant approach to solve this power control issue is to use a block code across the subcarriers 
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[13]. A well-designed code is able to provide a certain level of error protection and ensures a 
maximum PMEPR that is substantially reduced compared to uncoded transmission [12]. 

It has been proposed in a number of contributions [37], [36], [18] to use so-called Golay 
sequences [6] as codewords in OFDM since their PMEPR is at most 2 in g-ary phase-shift 
keying (PSK) modulation [27]. Major progress has been made by Davis and Jedwab [4]; slightly 
generalized by Paterson [23], it was proved that a large family of polyphase Golay sequences 
of length 2™ organizes in m!/2 cosets of a naturally generalized first-order Reed-MuUer code 
RMg(l, m) inside a generalized second-order Reed-MuUer code ZRMq(2, m) (see Definition [5]). 
A so-called supercode was then constructed by taking the union of these m!/2 cosets. However, 
due to a vanishing rate, these codes are only suitable for, say, m < 5. For this reason, Davis 
and Jedwab [4] proposed to include further cosets of RMg(l,m) with PMEPR greater than 2 
(but still low), in order to increase the rate of the codes at the cost of a slightly larger PMEPR. 
An exhaustive search technique was employed to identify such cosets. However, due to high 
complexity, such a search becomes prohibitive when m > 4 and q increases. 

Paterson's work [23] and the more general study by the author and Finger [32] provided 
some advanced theoretical background: it was shown that each coset of RMg(l,m) is entirely 
comprised of sequences lying in so-called complementary sets of the same size. It follows that 
the PMEPR of the codewords in each coset of RMg(l,m) is at most 2^^+^, where k depends on 
the algebraic structure of a coset representative. However this upper bound is not always tight, 
and at best it yields merely the lowest possible power of 2. 

Further results were obtained by Parker and Tellambura. In [21] they proposed another con- 
struction technique for cosets of RMg(l,m) comprised of sequences lying in complementary 
sets of small size. In [19], [20] the Rudin-Shapiro construction [33], [29] was used to obtain 
sequence sets with low PMEPR from a suitable pair of starting sequences (kernel). However such 
a construction generally does not yield cosets of RMq(l, m). In order to obtain complete cosets, 
an exhaustive search (although less complex than that proposed in [4]) for suitable kernels was 
performed in [19] and [20]. 

Still, some phenomena cannot be fully explained. For example Davis and Jedwab observed [4] 
that there are cosets of RM8(1, m) whose maximum PMEPR, taken over all words in the coset, 
is equal to 3. In addition numerical results suggest the existence of more coset classes whose 
PMEPR is bounded by low values not being a power of 2. This motivates further analyses of 
the PMEPR of cosets of RMg(l, m). In particular it is of interest to identify cosets with PMEPR 
close to 2. 

The remainder of this paper is structured as follows. In the next section we introduce the 
OFDM communication model and establish most of our notation. In Section |lll] we use the 
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Rudin-Shapiro construction to obtain sets of so-called near-complementary sequences from a 
given kernel. This construction is then stated in Theorem |7] in the context of algebraic normal 
forms of generalized Boolean functions. These preliminary results serve as a stepping stone to 
establish our main result, summarized in Theorem [10] and Corollary [151 where we construct 
cosets of RMq(l, m) from a kernel and prove an upper bound on the PMEPR of these cosets. In 
Section HVl we prove a new lower bound on the PMEPR of arbitrary cosets of RM^(l,m) and 
apply the results to those cosets constructed in Corollary [151 These observations lead to some 
relations between the PMEPR and the peak-to-average power ratio (PAPR) under the Walsh- 
Hadamard transform. In Section [V] we will comment on this issue and briefly discuss implications 
for the coding problem in multicode code-division multiple access (MC-CDMA) systems (cf. 
[24], [25]). 

The application of our theory is then illustrated in Section [VI] by providing some construction 
examples. First, it is shown that the Davis-Jedwab construction of Golay sequences [4] just 
arises as a special case in our theory, since they originate from trivial kernels (of length 1). 
Second, we present further classes of cosets of RMg(l, m) whose maximum PMEPR is between 
2 and 4. In particular we identify a class with PMEPR at most 3, and therefore, we provide a 
proof for the conjecture by Davis and Jedwab on the PMEPR of a subset of this class [4]. The 
upper bound on the PMEPR is also compared with our lower bound, and it is shown that the 
upper bound is attained in many cases. Finally, in Section IVII[ we close with some concluding 
remarks and open problems. 

II. Notation and Preliminaries 

Throughout this paper ^ = exp(V— 1 27r/g) denotes a primitive gth root of unity and q is an 
even positive integer. 

A. Aperiodic Correlations and Complementarity 

Let A,B e €^ be two sequences with A = {Aq Ai ■ ■ ■ An-i) and B = {Bq Bi ■ ■ ■ Bn-i). 
If Ai = and G Zg for alH = 0, 1, ■ ■ ■ , n — 1, we shall call A a polyphase sequence. The 
aperiodic cross-correlation of A and S at a displacement £ G Z is given by 

E Ai+.B* 0<£<n 

1=0 

C(A, B){1) = {"^Y, AB*_^ -n<£<0 
otherwise. 
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where ()* denotes complex conjugation. The aperiodic auto-correlation of A at a displacement 
£ e Z is then conveniently written as 

A{A){l)^C{A,Am. 
Definition 1: For two sequences A, B e C" we define the operator as follows 

n— 1 

A^B^Y. \A{Am+A{Bm. 

e=i-n 

If A and B are polyphase sequences of length n, we obtain 

n— 1 

e=i 

A pair of polyphase sequences {A, B) is called a complementary pair if A-k B — 2n, which 
implies that the aperiodic auto-correlations of the two sequences sum up to zero for each nonzero 
shift. In tribute to Golay, who extensively studied binary complementary pairs in connection with 
multislit spectrometry [6], they are often called Golay complementary pairs, and each sequence 
lying in such a pair is termed a Golay sequence. 

If A and B are polyphase sequences. A* B is at most 2n^. We shall call a pair of polyphase 
sequences {A, B) a near-complementary pair and the sequences therein near-complementary 
sequences if 2n < (A*B) <^ 2n^. In other words, the aperiodic auto-correlations of A and B 
sum up to small values at a few nonzero shifts. 

B. OFDM Power Control 

Let us consider an n-subcarrier OFDM system. The transmitted OFDM signal is the real part 
of the complex envelope, which can be written as 

n— 1 

S{A) {9) = 5^ Ai e v^2.(*+C)e^ < ^ < 1, 

i=0 

where C is a positive constant. The sequence A = {AqAi ■ ■ ■ An^i) is called the modulating 
codeword of the OFDM symbol. Throughout this paper it is assumed that the elements of A 
are selected from a PSK constellation. Then A is a polyphase sequence. 
The PMEPR of the OFDM signal (or of the modulating codeword A) is then defined to be 

PMEPR(A) ^ - sup \S{A){e)\^. (1) 
n o<e<i 

Notice that the PMEPR can be as large as n, which occurs, for example, if A is the all-one 
word. However it is desirable to use codewords with PMEPR that is substantially lower than n. 
For the construction of such codewords the following theorem will be essential in the sequel. 
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Theorem 2: Let (A, B) be a pair of polyphase sequences of length n. Then the PMEPR of 
A and B is at most {Ai^ B)/n. 

Proof: It is well known (cf., e.g., [35], [4]) that 

n-1 



\S{Am\' = A{Am + 2J2^ {AiAm e 




where 3?{.} is the real part of a complex number. Hence 

n-1 



\S{A){d)\' + \S{Bm\^ = 2n + 2 5^ ^{\A{Am + e 




£=1 
n-1 

< 2n + 2^|A(A)(£) + = A*S, 

^=1 

by Definition [U Since |S'(.)(6')p of each individual sequence is non-negative, the PMEPR of A 
and B is at most [^A-k B)jn. □ 
The above theorem is consistent with the results in [27] in the special case where (A, B) is 
a Golay complementary pair. Then the PMEPR of A and B is at most 2. In [37], [36], [18], [4] 
it has been proposed to exclusively use Golay sequences as codewords in OFDM. Consequently 
tight power control for OFDM is ensured, however, the code rate rapidly decreases for larger 
lengths. Theorem [2] motivates the use of larger sequence families with slightly higher PMEPR 
as codewords in OFDM. 

C. Generalized Boolean Functions and Associated Sequences 

A generalized Boolean function f is defined as a mapping / : Z™ Zg. Such a function can 
be written uniquely in its algebraic normal form, i.e., / is the sum of 2*" weighted monomials 

2^-1 m-1 

f{x) = f{Xo,Xi,--- ,Xm-l) = X] JJ X'^, 

where the weights cq, ■ ■ ■ ,C2™-i are in Zg, and (ioii ■ ■ -im-i) is the binary representation of 
< z < 2™, such that i = ^^Sq ij2^ is its binary expansion. The order of the ith monomial is 
defined to be J2^=o ij, and the order, or algebraic degree, of a generalized Boolean function f, 
denoted by deg(/), is equal to the highest order of the monomials with a nonzero coefficient in 
the algebraic normal form of /. 

A generalized Boolean function may be equally represented by sequences of length 2"^. We 
shall define the sequence 

^(/) = (/o/l---/2'"-l) 
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as the 'Lq-valued sequence associated with f and the sequence 

as the polyphase sequence associated with f. Here we denote = /(io,^i, • • • ,im-i)^ where 
(io^i ■ ■ -im-i) is the binary representation of < i < 2"\ 

In the remainder of this subsection we introduce the technique of extending polyphase se- 
quences of length 2™ and their corresponding generalized Boolean functions. It will be used in 
the next section to prove our results on near-complementary sequences. 

Definition 3: Suppose m > k, let f : ^ be a generalized Boolean function in the 
variables Xq, Xi, - ■ ■ , x^-i, and write F = ^{f). Let 

< io < ii < • • • < ik-i < TTi 

and write 

< Jo < ii < • • • < jm-k-i < m 

for the remaining indices. Also denote x = {xj^ ■ ■ ■ Xj^_^_^) and let d= [d^di - ■ ■ dm-k-i) be a 
binary word of length m — k. We define the extended sequence F\x=d\ of length 2"* as follows. 
As {uqUi ■ ■ -Uk-i) ranges over Z|, at position 

fc— 1 m— fc— 1 

a=0 a=0 

the sequence F[x=di is equal to ^/("0'"i'' '"^-i) and equal to zero otherwise. We also define the 
extended generalized Boolean function f\x\ that is formally regarded as a generalized Boolean 
function in m variables, i.e., it is of type — > Z^. This function is obtained from / by replacing 
the variables Xa by Xi^ in the algebraic normal form of / for a = 0, 1, • • • ,k — 1. 

Notice that F\^x=d\ comprises 2^ nonzero elements and 2"* — 2^ zeros. It is a consequence of 
the above definition that at the positions where F]^x=d] is nonzero the sequence F]^x=d\ is equal 
to the polyphase sequence associated with /[a,]. 

Example 4: Let / : Z| — > Z2 be given by 

Writing '+' for '+1' and '-' for we have F = *(/) = (+H h). Now take x = (xqXs) 

and <Z = (10). Then F^^=d\ = (0+0+00000-0+0000). By relabeling the variable indices in / 
according to 1— > 1, 1 1— > 3, we obtain the extended function f^^] — x^x^ + x^. Regarding this 
function as a generalized Boolean function in m variables, we obtain its associated polyphase 

sequence (+ + + + + + + -I \--\ 1-+). This sequence is equal to F\x=d\ at the positions 

where F\x=d\ is nonzero. 
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D. Generalized Reed-Muller Codes 

A code C of length n over the ring is defined as a subset C C Z^. Such a code is called 
Zq-linear if each Zg -linear combination of the codewords of C yields again a codeword of C. If 
C is Zy-linear, a co^e? of C is defined to be 

f + C^{f + c\ceC}, 

where / G Z^ is called its coset representative, and of course, the additions are taken modulo 
q. We say that a coset of a code Ci lies inside a code C2 if / + Ci C C2. 

We are interested in codes defined by generalized Boolean functions. In what follows we recall 
the definitions and some basic properties of the generalized Reed-Muller codes RM^ (r,m) and 
ZRMg(r,m) (cf. [4] and [23]). 

Definition 5: (a) For < r < m the code RMg(r, m) is defined as the set of sequences ^/'(/), 
where / is a generalized Boolean function Z™ Z^ of order at most r. 

(b) For g > 4 and 1 < r < m the code ZRMg(r, m) is defined as the set of sequences ip{f), 
where / is a generalized Boolean function Z^ Z^ with algebraic normal form containing 
monomials of order at most r — 1 and two times the monomials of order r. 

The codes RM^ (r,m) and ZRMg (r,m) are Zg-linear, and their minimum Lee distances are 
equal to 2™^'' and 2'""'"+^, respectively [4], [23]. In this paper we will particularly study cosets 
of the code RMg(l,m), which is comprised of the codewords corresponding to all affine forms 
over Zg in m two-state variables. Hence the number of words in such cosets is equal to 

III. Constructions of Sequence Families with Low PMEPR 
A. Rudin-Shapiro Constructions 

In what follows we present a slightly generalized version of the Rudin-Shapiro construction 
[33], [29], and exhibit its application to the construction of near-complementary pairs. Our main 
argument is the following lemma. 

Lemma 6: Let A and B be two sequences of the same length and let C = A + B and 
D = A-B. Then C D = 2(A^B). 
Proof: It is straightforward to show 

/1(C) (£) = A{A){(i) + A{B){e) + C{A, B){£) + C{B, A){C) 

A{Dm = A{A){1) + A{Bm - C{A, B){£) - C{B, A)(£). 

Combining the relations above and Definition [H the lemma follows. □ 
It is well known that the Rudin-Shapiro construction can be employed to recursively construct 
longer complementary pairs starting from a known complementary pair [6], [2]. Using the 
notation in the present paper, this will be illustrated in the context of a more general framework. 
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Let A and B be two polyphase sequences of length 2^. The extended sequence A[a.j.=o] is a 
sequence of length 2^+^ and contains the sequence A in the left half, while its right half contains 
only zeros. Likewise the sequence B[xi,=i] contains zeros in the left half and the sequence B in 
the right half. Observe that A(A[^,=o])(^) = A{A){i) and = A{B){e) for each 

i E li. Let us construct 

and notice that the same pair could be constructed by 

C={A\B), D = {A\-B), 

where (.|.) denotes concatenation. The reader may recognize the classical Rudin-Shapiro con- 
struction. By Lemma[6]we know that d^D = 2(Ai^B). Hence, if (A, B) is a complementary 
pair, so will be (C, D). This construction is also known as Golay's concatenation technique for 
synthesizing complementary pairs [6]. 
We may also find a new sequence pair by 

C = Ai^g=o] + Bi^g=i] 

D = A[^o=o] - -B[xo=i]- 

Notice that A[2.q=o] is obtained from A by inserting zeros in A at alternating positions starting at 
the second position. Likewise B[^o=i] is obtained from B by inserting zeros in B at alternating 
positions starting at the first position. It is then easy to verify that A[^o=o] ^B[^g=i] = Ai^ B, 
and by Lemma |6l we have C D = 2(A B). Therefore, if (A, B) is a complementary pair, 
so will be (C,D). This construction is essentially the same as Golay's interleaving technique 
for synthesizing complementary pairs [6]. 

If A and B are extended in more than one variable and Lemma [6] is applied multiple times, 
then we can obtain even more general methods to construct longer sequence pairs from short 
ones. Again, if we restrict our attention to complementary pairs, we obtain Golay's general 
interleaving method to construct complementary pairs of length 2*" [7], [6], [2]. 

We have now established the following. Starting from an initial polyphase sequence pair 
(A, B) of length 2'^ we can use the notion of extended sequences and a generalized Rudin- 
Shapiro construction to compose polyphase sequence pairs (C, D) of length 2"* with m > k. 
Since by Theorem |2] we have {C'kD)/2'^ = {Ai<B)/2^, the sequences A, B, C, and D have 
the same PMEPR upper bound (A -k B)/2'^. If the sequence pair (A, B) cannot be obtained 
from a shorter sequence pair in this way, then it is referred to as a kernel. 
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B. Explicit Constructions of Near-Complementary Sequences 

Using the language of generalized Boolean functions, an explicit construction for near-complementary 
sequences is presented in the following. 

Theorem 7: Let m > k and write m — k — s + t for non-negative integers s,t. Define the sets 

J — {0, - ■ ■ , s — 1,171 — t, ■ ■ ■ , m — 1} 

and / = 1tn\J . Denote the elements of J and / by 

< jo < ji < • • • < im-fc-i < -m 



and 



< io < ii < • • • < ife-i < m, 



respectively. Let a, 6 : ^ be two generalized Boolean functions and define / : — > Z^ 
by 

fi^Oj ■ ■ ■ T^m-i) = (i{x^^, ■ ■ ■ — Xj^^^^^) + b[xif^, ■ ■ ■ , Xii^_-^)xj^^g^ 

m—k—2 m—k—l 
a=0 a=0 

where Wq, • • • , Wm-k-i: w eZg and tt is a permutation of {0, 1, • • • , m — /c — 1}. Then 
and 

PMEPR(*(/)) < HW^. 

Proo/- Write cc = {xj^^^^ ■ ■ ■ Xj^^^_^_^^) and let d = (do • • • c^m-ifc-i) be a binary word of 
length m — k. Moreover let A — ^{a) and B — ^{b). Consider the extended sequences A[x=d] 
and B[x=d]- These sequences can be constructed by successively inserting zeros at alternating 
positions, at the beginning, or at the end of the sequences A and B. It is then straightforward 
to establish that 

A[x=di * B[x=di = Air B. (2) 

We now use a method similar to Golay's general interleaving construction, which is here 
applied to pairs that are not necessarily complementary. Let us define the recurrence formulae 

C^'^+^Hd^+l, • • • , dm-k-l) = C^''\0, d^+i, • • • , dm-k-l) + C'-D^'^Hh • • • , dm-k-l) 0) 

D^^'+^\d,+i, dm-k-l) = C(^)(0, ■ ■ • , dm-k-l) - C'D^^\l, d^+i, ■ ■ ■ , dm-k-l), 

(4) 



10 



where = 0,1, ■ ■ ■ ,m — k — 1 and 

C^^\dQ, di, - ■ ■ , dm-k-i) = Aix=d] 
D^^\do, di, - ■ ■ , dm-k-i) = B^x=d]- 

Since the positions of the nonzero components in C^^\0,d ^j^i, ■■■ ,dm~k-i) and 
D^^\l, d^+i, ■■■ , dm^k-i) are distinct for each /i = 0, ■ ■ ■ , m - A; - 1, C^™"^') and D^"^-^^ 
are polyphase sequences of length T^. Observe that y4(^"'''i:)(^))(£) = A{D''^''^){tj for each 
£ G Z. Therefore, by repeated application of Lemma [6] and using Q, we have 

^{m-fc)^^{m-fe) _ 2'"-'^(A*B). (5) 

We know that A\^x=d] and the polyphase sequence corresponding to a\^x] are equal in those 
positions where A\^x=d\ is nonzero. An analogous statement holds for B[x=d] and h[x]- Similarly 
we can find functions, say c*^'^^ and d^^\ whose associated polyphase sequences are equal to C^^-* 
and I?*-^'^ at their respective nonzero components. Using ([3]) and (|4]) it can be verified that the 
functions c^^'> and rf*^^^ can be recursively constructed as follows 

where c^°) = a[x] = o-ixi^, ■ ■ ■ , Xt^^^J and = 6(3;] = b{xig, ■ ■ ■ , Xi^_^). Explicitly we obtain 

.(/^) _ 

fi—2 fJ-—^ 



a=0 a=Q 



where /x = 1, 2, ■ • • ,m — k. We have / = c*^"^"^) + w and / + f = d'^^~^^ + w. The 
theorem follows then from ([5]) and Theorem [21 □ 
We refer to Section IVIII for a discussion on the relation of Theorem |7] to previous results in 
[19] and [20]. 



C. Cosets o/RMg(l,m) with Low PMEPR 

Once suitable kernels are known. Theorem |7] identifies a large family of sequences with 
low PMEPR. However it would be desirable to construct sequence families that naturally form 
unions of cosets of RMq(l, m) inside a higher-order generalized Reed-MuUer code. In this way 
we could quickly obtain error-correcting codes, for which well-designed encoding and decoding 
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algorithms exist (see, e.g., [4], [26], [8], [30], [31]). The remainder of this section is dedicated 
to finding such sequence sets. First we require some preliminaries. 

Definition 8: Let / : — * be a generalized Boolean function. We define the sequence 
$(/) of length (4^ + 2)/3 as follows. As (uqUi ■ ■ ■ Uk-i) ranges over Z2, at position 

k-l 

the sequence $(/) is equal to ^/("o,«ir - ^nd equal to zero otherwise. 

We remark that the sequence $(/) may also be obtained from the extended sequence '^{f)[x=d\ 
by setting x = (xi X3 ■ ■ ■ X2k-3), letting d be the all-zero word of length k — 1, and deleting the 
trailing zeros. 

Definition [8] implies that for any generalized Boolean function / : Z2 — Z^ at its 2^^ nonzero 
elements the sequence $(/) coincides with the polyphase sequence associated with a function 
that is obtained by replacing Xa by X2a for each a = 0,1,-- - ,k — 1 in the algebraic normal 
form of /, where this new function is regarded as the algebraic normal form of a generalized 
Boolean function in 2 A; — 1 variables. 

Example 9: Let / : Zg ^ Z4 be given by 

/(xo, Xi) = 2xqXi + 3xq + Xi. 

We obtain: 

$(/) = (! -J OOj -1), j = V^. 

At the nonzero positions the above sequence coincides with the polyphase sequence associated 
with the function (7 : Z2 — Z4 whose algebraic normal form is given by g{xo, xi, X2) = 2xoX2 + 

3Xo + X2. 

We are now in the position to state the main theorem of this paper. 

Theorem 10: Suppose that m > k. Let a, 6 : Z^ ^ Z^ be two generalized Boolean functions, 
and let / : Z™ Z^, a generalized Boolean function in the variables xq, Xi, ■ ■ ■ , x^^i, be given 
by 

/(Xo, ■ ■ ■ , Xm-l) = a(a;,r(0), " " " ? a^7r(fc-l))(l — a^7r(fc)) + K^niO), ■ ■ ■ , XTr{k-l))Xn{k) 

+ I ^ Xn{a)XT,{a+l) + ^ WaXn{a) + 
a=k 0=0 

where tt is a permutation of {0, 1, ■ ■ ■ , m — 1} and wq, ■ ■ ■ , Wm-i,w E Zg. Then 
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and 

PMEPR(^(/)) < ^ '^^ ^ ^ 

We need a series of lemmas in order to prove the theorem. It is well known that each integer 
< i < 2"^ has a unique binary representation {iq - • •im.-i)'> such that i — Y^^Zq ia'^°' 
la e {0, 1}, If we now allow e {—1, 0, +1}, we obtain a signed-digit representation (SDR) 
of i [11]. Such a representation is in general not unique. An SDR is called sparse if it contains 
no adjacent nonzero entries. We have the following lemma, which is a well-known result in 
number-representation theory (cf. [9], [11]). A proof is included, since the lemma will play a 
crucial role in the sequel. 

Lemma 11: Every nonzero integer i has a unique sparse SDR {iq - ■ ■ im-i) with im-i 7^ 0. 
Proof: Without loss of generality we assume that i > 0. We first show that for each nonzero 
integer there exists a sparse SDR. For i — 1 this is obvious. The remaining cases are proved by 
induction. There are three cases: (i) if i is even, take a sparse SDR for i/2 and prepend '0', (ii) 
if i = 1 (mod 4), take a sparse SDR for (i — l)/4 and prepend '1, 0', (iii) if i = — 1 (mod 4), 
take a sparse SDR for and prepend '—1,0'. 

Now let us prove the uniqueness of a sparse SDR for a given positive integer. We first show 
that there is a unique sparse SDR for i — 1 and then proceed by induction. Suppose im-i — 1- 
Then 



i > < 



2"^-i - y 2^" = — — if m is odd 
^ 3 

a=0 

2^ + 2 

2m-i _ J2 22"+i = — — if m is even. 



a=0 

If ijn-1 = — 1, a similar argument yields i < 0. We conclude that, if i — 1, then m — 1 and 
iq — 1. Hence a sparse SDR for i = 1 is unique. 

Now suppose that i > 1 is the smallest positive integer with two sparse SDRs, namely 
(wo, ■ ■ ■ , Ujn-i) and {vq, . . . , Vn-i). If i is even, then uq — Vq — 0. By shifting the two SDRs one 
position to the left, we obtain two sparse SDRs for i/2 < i, which is a contradiction. If i = 1 
(mod 4) or i = —1 (mod 4), then uq — vq — 1 or uq — vq — —1, respectively. Consequently, 
we would have two sparse SDRs for (i — uo)/4 < i. Again we arrive at a contradiction, which 
completes the proof. □ 

We are now able to prove the following three lemmas. 

Lemma 12: Suppose a, 6 : Z| — > are two generalized Boolean functions. Then 
holds for any appropriate list of variables x and for any d of suitable length. 
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Proof: Let 

< zo < "^1 < • ■ • < ^fc-i < ^ 

and write 

< jo < ji < ■ • • < jm-k-i < m 

for the remaining indices. Denote ) and let d = {cIq - ■ ■ dm-k-i) be a fixed 

binary word of length m-k. Write A = '^{a)^^^^], A = ^{a), B = '^{b)[x=d], and B = ^{b). 
We claim that 

AiAm = J2 MAW) (6) 

i'ec{i) 

for some disjoint sets C{i), which we shall now prove. Let u and u' have binary expansion 

k—l m—k—l 

«=0 a=0 
k-l 

•,2a 



U 

respectively, and use an analogous notation for v and v' . Then u and u' define the positions of 
the nonzero components of A and A, respectively. By Definition [3] and Definition [8] we have 
Au = Au' for any u and corresponding u'. Now consider the nonzero product Au'Al, = A^Al 
with u' 7^ v' (and therefore u ^ v) occuring in the expression A{A){u' — v') and also in 
A{A){u - v). An SDR of d' = u' - v' is given by 

{uq-Vo Ui-Vi---0 Uk-l-Vk-l). 

This SDR is sparse and, by Lemma [TTl unique. In other words, for each a = 0, ■ ■ ■ , — 1 the 
differences Ua — Va are uniquely determined by £'. Hence £ = u — v (which is independent of d) 
is also uniquely determined by This means that all nonzero products contributing to the sum 
in the expression A(A){i') also contribute to the sum in the expression A(A){£) for exactly one 
i, which proves Q. 

Now we use Q to establish 

A^B = 2^+^ + 2 J2 + A{B){e)\ 



1=1 

2^-1 



)fc+l 



A(A)ir)+ J2 MB){e] 



1=1 

2'"-l 

<2'=+i + 2^ \A{A){e) + A{B){e 

1=1 
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Since the sets are disjoint and each product contributing to A(A)(i') for some £' also 
contributes to A{A)(i) for some i, the latter expression is equal to Ai<B, which completes the 
proof. □ 
Lemma 13: Let a, 6 : Z2 — be two generalized Boolean functions, and define a = a + L 
and b = b + L, where 

fe-i 

L = ^WaXa, Wo, ■ ■ ■ , e Zg. 

a=0 

Then 

$(a)*<l>(6) = <l>(a)*$(6). 

Proof: Write A = $(a), A = $(5), B = $(6), and B = $(6). Let u and w have binary 
representation {uq ■ ■ ■ U2k-2) and [vq ■ ■ ■ V2k-2), respectively. For u ^ v we consider the products 
appearing in A{A)(u — v) 

A^A* =Au^^"^°^"^'^" A*^~'^"=o'^"^'^" 



Let us focus our attention on the word (mq — Vq - ■ ■ U2k-2 — ^2^-2) » which is an SDR of u — v. 
By Definition [H Au is equal to zero if and only if there exists an a E {0, ■ ■ ■ , k — 2} such 
that U2a+i 7^ 0. Therefore, if the product AuAl is nonzero, an SDR of m — is sparse and, 
by Lemma [TTl unique. In this case the differences U2a — V2a are uniquely determined by the 
difference u — v. Then we have A^A* = A^A*^-^*^""^), where K{u — v) only depends on m — t> 
and not explicitly on u and v themselves. Thus 



and similarly. 



We conclude that 



A(A)(£)=e^WA(A)(£), 



\A{Am+A{Bm\ = \A{Am+A{Bm\ 



for all i E and the lemma follows. □ 
Lemma 14: Let a, 6 : Z2 ^ be two generalized Boolean functions and define 

a(xo, ■ ■ ■ ,Xk-i) = a{xa(o), ■ ■ ■ ,Xa(k~i)) 
b{xo, ■ ■ ■ = 6(x^(o), ■ ■ ■ 
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where a is a permutation of {0, 1, ■ ■ ■ ,k — 1}. Then 

Proof: Write A = $(a), A = $(a), B = ^{h), and B = ^(b). Let u and m' have binary 
representation 

{UqUiU2 - ■ ■ U2k-ZU2k-2) 

and 

(m2(t(0) «1 U2a(l) ■ ■ ■ U2k-3 «2<T(fe-l)), 

respectively. In an analogous manner we denote the binary representations of v and v'. Let us 
consider the product AuAl = Au^A^,. By the same reasoning as in the proof of Lemma [13] we 
know that, if Au'Al, is nonzero, an SDR of u' — v' is sparse and uniquely determined by u' — v'. 
Similarly an SDR of m — is sparse if is nonzero. It is then clear that, if Au'A*^, is nonzero, 
there is a one-to-one correspondence between the sparse SDR of u — v and the sparse SDR of 
u' — v' , and hence, u — vh uniquely determined by the difference u' — v' and does not explicitly 
depend on u' and v' themselves. This means that each nonzero product appearing in A{A){u' —v') 
also appears in A{A){u — v) and, by setting a := o^^ and arguing analogously, the converse 
holds too. Hence A(A){u -v) = A{A){v! - v') and, similarly, A{B){u - v) = A{B){u' - v'). 
The lemma follows then immediately. □ 
Proof of Theorem Ud\- Let a set of m — A; indices be given by {jo, Ji, " " " , jm-fe-i}: where 
< ja < for each < a < m — k. Write x = (xj^ ■ ■ ■Xj^ ^. J and let d be a binary word 
of length m — k. First we construct a number of initial function pairs (a, b), where 

a{xo, ■ ■ ■ , Xk^i) = a(x<^(o), ■ ■ ■ , x^(^k~i)) + L{xq, ■ ■ ■ , Xk-i) 
b{xo, ■ ■ ■ = &(x^(o), ■ ■ ■ ,x^(^k-i)) + L{xo, ■ ■ ■ ,Xk-i), 

cr is a permutation of {0, 1, ■ ■ ■ ,k — 1}, and 

fc-i 

L{Xo, ■ ■ ■ , Xk-l) = ^ WaXa, Wq, ■ ■ ■ , Wk-1 E Zg. 

From Lemma [121 Lemma [131 and Lemma [14| it follows 

^(a)[.=d] * ^(&)[.=d] < Ha) ^ $(6). (8) 

We can then proceed with the same constructive reasoning as in the proof of Theorem [7] to 
obtain the function pair (c^"^~''\ from (a, b). Using ([8]) and the arguments from the proof 

of Theorem [7] we conclude 

^ m{d^""^^) < 2'^~'=[$(a) ^ $(6)]. (9) 
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Then c*^™^''^ is a function where the addition of all linear terms in m variables is possible. 
Moreover any permutation can be applied to the indices of the m variables in 0*^"*"*^^ because a 
and the indices jo,ji, ■ ■ ■ ,jm-k-i can be chosen arbitrarily, so that all < Jq, < m are distinct. 
We can, therefore, replace ja by 7r(k + a) for each a = 0, 1, • ■ ■ , m — k — 1 and set / = c^"^~'''^+w 
and / + = d'^"^^^^ + w. The theorem follows then from ^ and Theorem |2l □ 

We have a useful corollary of Theorem [lOl 

Corollary 15: Let m > k and suppose a, 6 : Z| ^ Zq are two generalized Boolean functions 
in k variables. Then each polyphase codeword in the cosets ^/'(/) + RMq(l, m) with / : Z^ 
given by 

m—2 

f = fi^O, ■ ■ ■ 5 Xm-l) = - a^7r(a)3;7r(a+l) 

a=k 

+ a.(x^(o), ■ ■ ■ , a;^(fc_i))(l — a;,r(fc)) + K^wio), " " " , a;7r(fe„i))a;7r(fe), 

where vr is a permutation of {0, 1, ■ ■ ■ ,m — 1}, has PMEPR at most [$(a) $(6)]/2^ and lies 
inside RM^ (r,m) with 




max{deg(6 — a) + 1, deg(a)} if m = k + 1 
max{deg(6 — a) + 1, deg(a), 2} if m > k + 1. 



In particular, if g > 4 and all coefficients of the monomials in the algebraic normal form of 
/ with degree equal to r are even, then the cosets are contained in ZRMg (r,m). The number 
of distinct functions /, and therefore the number of distinct cosets, is at least (m — k)\/2 and 
at most m! (the true value depends on how many permutations can be applied to the variable 
indices in the algebraic normal forms of a and b such that distinct pairs (a, b) are generated). 

In order to construct families of cosets of RMg(l, m) with low PMEPR, we just have to find 
two generalized Boolean functions a, 6 : Z2 — > Zg with [$(a) <l>(6)]/2'^ being low. This can be 
accomplished by an exhaustive search and, if k is small enough, such functions could even be 
found using a 'by hand' construction technique. Several examples will be given in Section |Vll 

IV. Lower Bounds on the PMEPR 

Several lower bounds on the maximum PMEPR taken over all the words in a coset of 
RMg(l,m) have been proposed in [23], [34], [16]. These approaches rely on the examination 
of the OFDM signal at time 6^ = or at some other 9. Initially this method was proposed 
in [3]. However existing results apply to second-order cosets of RMg(l,m), where the coset 
representative is binary (in the g-ary context this means that it has values '0' and 'g/2'). 
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In the following we provide a general lower bound on the PMEPR of the cosets ip{f) + 
RMq(l,m), where / : Z™ Zg is an arbitrary generalized Boolean function. A g-ary general- 
ization of the Walsh-Hadamard transform plays an important role in our study. 

Definition 16: We define the g-ary Walsh-Hadamard transform (WHT) of a generalized Boolean 
function / : Z™ Z^ (or of its polyphase sequence to be 

F{w) = 

where w E Z™ and '■' denotes the scalar product of vectors. By convention, in the special case 
where m = and / = a (a G Zg), we define F to be equal to 

If g = 2, the above definition coincides with that of the classical WHT (see, e.g., [15, Chapter 
14]). The following lemma will be useful in the sequel. 

Lemma 1 7: Let /, g : Z^ Z^ be two generalized Boolean functions that are related by 

g(x) = f{Ax) + V ■ X + V 

with V E Z^, V E Zg, and A being an mxm permutation matrix. Then the sets {G{w) \ w E Z™} 
and {^'"F{w) \ w E Z™} are equal. 
Proof: Write 



;f(Ax)+w-x+v-x+v 
xeZ^ 

Setting X = A~^y, we obtain 



= CF{A{w + v)), 

and the lemma follows. □ 
Remark: In the particular case where q = 2, Lemma [IT] still holds for A being an arbitrary 

invertible binary matrix (cf. [15, Chapter 14]). 

Now we are ready to formulate the following theorem. 

Theorem 18: Let / : Z™ ^ Zg be a generalized Boolean function. Then there exists a 
polyphase codeword in the coset ip{f) + RMg(l,m) having PMEPR at least 

— max \F(w)\^. (10) 
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Proof: Consider the coset + 'Ryiq{l,m), which contains the polyphase codewords 
Fww = ^ifix) + w-x + w), where w E and w E TLq. The complex envelope corresponding 
to -F,,„„ reads for ^ = 



WW 

2'"-l 



>S'(-F'i£;^) (0) ^ ^ PwW, 

1=0 



i=0 

f(x)+w-x+w 



Therefore 



max sup |^(F^J(e)|' >max|^(F^J(0)|' 

'^'^ 0<9<1 

= max|F(io) 

w,w 

= max \F{w)\'^, 

w 

and the theorem follows from ([T]). □ 
Remark: In the case where / is a quadratic Boolean function and g = 2, it can be shown that 

[16, Theorem 1] coincides with the above theorem. 

General lower bounds on the expression in (fTOl) directly give lower bounds on the achievable 

PMEPR of cosets of RMg(l, m). Let us first discuss the binary case, i.e., g = 2. It is well known 

(cf., e.g., [15, Chapter 14]) that 

max|FM|=2--2 min d^, c) , 

idGZ™ ceRM2(l,m) 

where d//(-, ■) is the Hamming distance between two binary sequences. The expression 

pirn) = max min dn (/, c) 

f&l"^ ceRM2(l,m) 

is the covering radius of RM2(l,m). Therefore, provided that p(m) < 2™~\ we have 

min max = 2™ — 2p(m), 

where the minimum is taken over all Boolean functions / : Z2. Results on the covering 

radius of RM2(l,m) can now be used to lower-bound (flOl) for q = 2 (cf., e.g., [1], [25]), 
and therefore, to lower bound the PMEPR of cosets of RM2(l,m). If m is even, we have 
p(m) = 2™^^ — 2™/^^^, which leads to the trivial lower bound of 1 for (fTOl) . For odd m it is 
known that p(m) > 2™^^ — 2*^™^^^/^, where equality holds when m < 7. We conclude that the 
PMEPR of cosets of RM2(1, m) is at least 2 if m is odd and m < 7. For odd m > 9, the exact 
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value of p(m) is unknown. However, if m > 15, it is known that p(m) > 2™^^ — ||2'^™~^)/^, 
which implies that (flOl) can be strictly smaller than 2 in this case. 

Much less can be said when q > 2. However at least we can assert that the PMEPR of cosets 
of RMg(l, m) must be strictly greater than 1 if g is a multiple of 4. To this end observe 

max \F(w)\ > max \F(w)\. (11) 

It can be shown that the set {F(w) \ w E j'^a} is equal to the set of values of the 2"* possible 
{H, N}"^ transforms [28] of a function / : Zg. In [28] it was proved that there is no 

function / : — > Zg for which all 2"* {H, N}"^ transforms contain only absolute values 2™/^. 
Moreover we conjecture that at least for small values of m and for q being a multiple of 4 the 
right-hand side in (fTT]) is lower-bounded by 2'^"'+^^/^, and therefore the PMEPR of cosets of 
RMg(l,m) is at least 2. 

In what follows we apply Theorem [T8] to obtain lower bounds on the PMEPR of the cosets 
of RMq(l,m) constructed in the previous section. 

Theorem 19: Let a, 6 : Z^ ^ Zg be two generalized Boolean functions, and let / : Zg 
be given by 

q '""^ 

/'(•^O; ' ' ' ; -^m— l) ^ ^ •^7r(a)'^7r(a+l) 

a=k 

+ o(a;7r(o), ■ ■ ■ , a^7r(A;-i))(l — XT,(^k)) + &(a^7r(o), ■ ■ ■ 5 a;^(fc-i))a;7r(fc), 

where vr is a permutation of{0,l,---,m — 1}. Then there exists a polyphase codeword in the 
coset 'ijj{f) + RMg(l,m) having PMEPR at least 

^ max U(t(;')(l + 0+0 + -O+OT (^2) 



if m — A; is even and at least 

-1^ max \A{w') + B{w')r'\' (13) 
if m — A; is odd. Here A(w') and B{w') are the g-ary WHTs of a and 6, respectively, and 

W' = {wq ■ ■ ■Wfc_i). 

Proof: We intend to find a lower bound for 

1 



— max \F(w' 

2m ^yg^m I ^ ' 



|2 



where F(w) is the g-ary WHT of /. Using Lemma [17] we conclude that, in order to prove the 
theorem, it is sufficient to assume that vr is the identity permutation. So we are interested in the 
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coset representative corresponding to 

a=k 

Denote the g-ary WHT of fm by Fm{w). We require the following expansion 

where v = {wq ■ ■ ■ Wm-2)- Let us first consider the case m = /c + l, som — Zeis odd. We have 
fk+i = a{xo, ■ ■ ■ ,Xk-i){l - Xk) + fc(xo, • • • ,Xfc_i)zfc, 

and thus, 

/fc+i(xo, ■ ■ ■ , 0) = a(xo, ■ ■ ■ , Xk-i) 
fk+i{xo, ■ ■ ■ , Xk^i, 1) = b{xo, ■ ■ ■ , Xk-i). 

With (O we obtain 

Fk+i{w',Wk)=A{w') + B{w')r'. (15) 
Now let m > A; + 1. Then we have 

fm{XQ, ■ ■ ■ , Xm_2, 0) = fm-l{Xo, ■ ■ ■ , Xm-2) 

fm{Xo, ■ ■ ■ , Xm-2, 1) = fm-l{Xo, ■ ■ ■ , Xm-2) + ^2;m-2, 

and with (fUl) it follows that 



Fm{Wo, ■ ■ ■ ,Wm-l) = Fm-l{Wo, " " " , Wm-2) + -^^-1(1^0, " " " , Wm-2 + -^)C"" (16) 

For m — k odd suppose 

Fmiw', Wfe, 0, ■ ■ ■ , 0) = 2('"-'=-^)/2^AK) + Biw')^), (17) 

which is true for m = k + 1 (see (fTSi) ). We will use this expression as a hypothesis for the 
following induction. We employ (fT6l) to obtain 

F^+i(tf;', Wfc, 0, ■ ■ ■ , 0, w^) = 2(™-'=-i)/2 [A(t(;')(l + e""") + Biw')dl - T™)] (18) 

and 

F^+2(tf^', Wfc, 0, ■ ■ ■ , 0, w^, Wm+i) = 2(— [A(^') (1 + C"-) + Biw')^! - C"-) 

[A{w'){l - C"') + B{w')C''{l + r™))] • 
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Consequently we have 

F„+2(^', Wk,0,--- ,0) = 2(— {A{w') + B{w')C') 



which proves by induction that (frTI) and (|T8l) hold in general. Now we can write for m — A; even 

1 2 1 I ,2 

— max\Fraiw)\ >— max \Fra{w',Wk,0,--- ,0,Wk+i)\ 

^ max, I + 0+') + - 0+')!^ 
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and for m — k odd 



— max > — max iFmiw' ,wu,0, ■ ■ ■ ,0)^'^ 



^ max \A{w') + B{w')C"''? . 



Then the statements in the theorem follow from Theorem [181 □ 



V. The PAPR under the Walsh-Hadamard Transform 

Suppose that / : — > is a generalized Boolean function and p is a divisor of q. Then 
we define the peak-to-average power ratio (PAPR) of ^(/) under the p-ary WHT to be 

PAPRp(vl/(/)) A 1 n,ax \F{w)\\ 

2™- itiG(g/p)Z™ 

If p is omitted, we shall refer to the PAPR of ^(/) under the classical (2-ary) WHT. It is 
apparent that 

PAPRp(vI/(/)) < PAPR,(vI/(/)). 

We remark that this measure in fact arises from a more general definition of the PAPR under 
unitary transforms (cf. [21], [28]). Also the PMEPR can be restated in this context [21]. The 
PAPR under the p-ary WHT is an important measure in cryptographic applications [22], and it 
is of interest in MC-CDMA communications systems [25], particularly when p = 2. 

OFDM and MC-CDMA enjoy several similarities: in both cases codewords are used to 
modulate simultaneously a number of orthogonal signals, which are continuous in OFDM and 
discrete in MC-CDMA. In MC-CDMA systems the PAPR (under the classical WHT) turns out 
to be an analogous measure of the PMEPR in OFDM systems [24]. In [24] and [25] Paterson 
studied binary codes whose codewords have low PAPR, whereas constructing nonbinary codes 
is left as an open problem. 

Theorem [TSl relates the PMEPR of cosets of RMg(l,m) with the g-ary WHT of its coset 
representative. It is a consequence of Lemma [TTI that every word in a coset of RMq(l,m) has 
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the same PAPR under the g-ary WHT (but generally not under the p-ary WHT). Using this fact, 
Theorem [T8] states that cosets with low PMEPR must have also low PAPR under the g-ary WHT, 
and therefore, also under the p-ary WHT. Thus we have established the following corollary. 

Corollary 20: Let / : ^ Zg be a generalized Boolean function, and suppose that p 
divides q. Then the PAPR under the p-ary WHT of the polyphase codewords in the coset ip{f) + 
RMg(l, m) is at most the PMEPR of the polyphase codewords in the coset. 

The above corollary shows that all codes with low PMEPR that are unions of cosets of 
RMg(l, m) (e.g., codes that arise from the results in the present paper and the codes constructed 
in [4] and [23]) also enjoy low PAPR. Hence with Corollary [20] a number of (generally nonbinary) 
coding options for MC-CDMA with low PAPR can be derived. However, as opposed to many 
of the codes in [24] and [25] having PAPR equal to 1, Corollary [20] in combination with known 
results (e.g.. Corollary [T5]) yields an upper bound on the PAPR that is at least 2. 

We refer to the final section of [25] for further discussions on the connection between codes 
for OFDM and MC-CDMA. 

VL Construction Examples 

In the following we will apply Corollary [15] in order to find families of cosets of RMq(l, m) 
with low PMEPR. By Corollary [20] such cosets also have low PAPR under the p-ary WHT, where 
p is a divisor of q. Trivially we may choose our kernel functions a and b from the complete set 
of generalized Boolean functions in k variables. In such a way, we can construct large sequence 
families with their PMEPR bounded by 2'''+^, since the trivial bound states 

2k - 

However it can be shown that these sequence sets are just subsets of a larger family of comple- 
mentary sets of size 2^"'+^ as constructed in [32] and [23]. Hence this approach does not yield 
any new sequence sets. 

A. Sequence Sets with PMEPR at most 2 

An immediate consequence of Corollary [15] (and in this case also of Theorem [T]) is the well- 
known construction of Golay sequences over Z^, which has been stated by Davis and Jedwab 
in [4] for q being a power of 2. 

Corollary 21: Each of the m!/2 cosets of RMq(l, m) with coset representatives corresponding 

to 

m-2 

-^^x^^i^x^i+i), (19) 

i=0 
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where tt is a permutation of {0, 1, ■ ■ ■ ,m — 1}, is entirely comprised of Golay sequences. In 
particular the maximum PMEPR of the polyphase words in these cosets is (i) exactly 2 if m is 
odd and q arbitrary (even) and if m is even and g = (mod 4), (ii) at least 1 + cos^ ^ and at 
most 2 if m is even and g = 2 (mod 4). 

Proof: The corollary follows by setting A; = 0, a = 0, and 6 = in Corollary [15] and 
Theorem [T9l Since ^{a)-k^(h) = 2, the upper bound on the PMEPR is an immediate consequence 
of Corollary \T5\ Note that a and b have g-ary WHT equal to 1. For odd m the lower bound on 
the PMEPR follows then immediately from (fT3l) . For even m and q = (mod 4) put 

-g g^ 



-4' 4 



in the maximum calculation in (fT2)) . which leads to a lower bound of 2. For even m and g 
(mod 4) it is not hard to see that the maximum in (fT2l) is attained for 

'g-2 g + 2^ 



w 



Hence, with j = v^^, ([12)) becomes 



\ max 1(1 + + ^"(1 - = 7 1(1 + Je^''^'') + je"^'^/'?(l - 

4 iuez2 ' ' 4 ' ' 



1 
4 



2 ^ 

1 + cos — , 



2 + 2j3fJ{e^"/'?}| 



which completes the proof. □ 
We remark that the upper bound is well known: it was first proved by Davis and Jedwab for 
g being a power of 2 [4, Corollary 6, Corollary 9] and generalized to arbitrary g by Paterson 
[23, Corollary 11]. The lower bound in Corollary [2T] (i) has been first proved Cammarano and 
Walker [3], and the case when m is odd also by Paterson [23, Theorem 21]. Except for the 
trivial bound 1 when g = 2, the lower bound in Corollary [2T] (ii) is new. Note that these upper 
and lower bounds on the PMEPR now arise in a uniform way from a general framework. From 
the discussion following Theorem [18] we also conclude that the considered cosets are optimal at 
least for g = 2 and odd m < 7 in the sense that for these parameters there cannot exist cosets 
with PMEPR lower than 2. 

Li and Chu [14] reported 1024 additional quaternary (with elements in {1,-1, V — 1, — v^— T}) 
Golay sequences of length 16. While Fiedler and Jedwab [5] recently provided an explanation 
for these sequences, it has been observed earlier by Holzmann and Kharaghani [10] that there 
exists essentially one quaternary kernel of length 8 

((+ + + - + + -+), + t t-)), J = V^,i = -J. 
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These two sequences are associated with generalized Boolean functions Z4 corresponding 

to the forms 

2xoXi + 2xiX2 and 2x9X2 + 2xiX2 + xq + xi, 

respectively. Following [10] from this kernel 512 equivalent kernels can be generated by applying 
several operations that preserve the complementary property. These kernels can be used in 
conjunction with Theorem |7] to obtain an explicit construction for Golay sequences of length 2'", 
where m > 3. Notice that all of them correspond to cubic forms. Indeed, if m = 4, we obtain 
exactly the 1024 Golay sequences of length 16 reported by Li and Chu [14]. Unfortunately, 
in this way, we cannot construct complete cosets of RMg(l,m) whose PMEPR is bounded by 
2. Instead Corollary [15] reveals that the PMEPR is at most 5 if the 512 kernels are used in 
conjunction with Corollary [15] to build cosets of RMq(l,m). 

B. Sequence Sets with Maximum PMEPR between 2 and 4 
Consider the functions a, 6 : Z2 ^ with 

a = a(xo,xi) = |xoXi 
h = 6(xo, xi) = -XqXi + [a + -j Xq + pxi, 
where a,P E {q/p)Zp and p divides q. We have 

$(a) = (11001 - 1) 
$(6) = (1 - C^f^) 



and 



(/l($(a))(£) |0<^< 6) = (4 00 1 - 1) 
(A($(6))(£) |0<£< 6) = (4 00 - e'^-" 



Hence 



$(a)*$(6) = 2(4+ |1 -^/^-"l + |1 < 16. 



So, by Corollary [T5l the PMEPR of the cosets of RM5(1, m) (m > 2) with coset representatives 
corresponding to 

m-2 
i=0 

is at most 4. These cosets lie inside RMg(2,m) and in particular inside ZRMq(2,m) if q > p. 

There are p"^ choices for the pair a, (3 and for each pair we can construct m!/2 distinct cosets 
by applying to the variable indices m!/2 permutations that are invariant under reversal. However 
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setting a = [3 = Q ox a = (3 = q/2 in (|20|) will generate the same set of coset representatives. 
Hence the total number of distinct coset representatives amounts to (p^ — l)(m!/2). 

For particular values for a and (3 the upper bound on the PMEPR of the cosets can be tightened. 
Ifa = /3 = 0orQ; = /3 = g/2, we obtain sequences with PMEPR at most 2. However in this 
case the coset representatives are of the form given in (fT9l) . More generally, the {p^ — l)(m!/2) 
cosets of RMg(l,m) organize in classes with the same PMEPR upper bound. For p = 2 we 
obtain just two classes: m!/2 cosets with maximum PMEPR 2 and m! cosets with PMEPR at 
most 4. Using Theorem [T9] it can be shown that the latter bound is tight if g = (mod 4). 

For g > 4 and p = 4 we have a, /3 G {O, f , |, ^} and the classes are: 



13 + a 


[3 — a 


max. PMEPR 


# cosets 








2 


m!/2 


arbitrary (0) 


(arbitrary) 


3 


2m! 


! (mod 1) 


1 (mod 1) 


2 + V2 


4m! 


arbitrary 


arbitrary 


4 


m! 



Notice that the number of cosets are counted by considering only those cosets that are not 
contained in a class with a lower PMEPR bound. For example the total number of cosets with 
PMEPR at most 3 is equal to 5m!/2, while there are 13m!/2 cosets with PMEPR at most 
(2 + \/2). Employing Theorem [19] it can be shown that the PMEPR bounds given above are 
tight if g = (mod 8). If g = 4, Theorem [T9] yields a lower bound of 2 for the first class, 2.5 
for the second and third class, and 4 for the fourth class. 

We remark that the second class includes those 48 cosets of RM8(1,4) inside ZRM8(2,4) 
for which Davis and Jedwab observed that their PMEPR is exactly 3 [4, Tables V and VI]. 
While Nieswand and Wagner [17] could explain this behaviour partially, a complete proof arises 
as a special case of Corollary [15] and Theorem [T^lfl . Moreover we identified further cosets of 
RMg(l, m) with PMEPR at most 3. 

For g > 8 and p = 8 there are already 9 classes with different PMEPR upper bounds. These 
are: 

'The upper bound for this particular case has been proved alternatively by the author and Finger in [32]. 
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max. PMEPR 


# cosets 




2 


m!/2 


2 + 1/v^ ^ 


2.707 


4m! 


2 + V 2 - V2 ^ 


2.765 


4m! 




3 


2m! 


2 + -^1 + 1/72 ^ 


3.307 


8m! 


2 + V2 ^ 


3.414 


4m! 


3 + 1/V2 ^ 


3.707 


4m! 


2+^2 + 72 ^ 


3.848 


4m! 




4 


m! 



Applying Theorem [T9l we are able to show that the PMEPR bounds are tight if g = (mod 16). 

The number of suitable kernels of length 4 producing cosets of RMg(l,m) with maximum 
PMEPR between 2 and 4 can be increased further. To this end, consider the pairs (a, 6), where 
a,b : Zl ^ Zg are given by 

a = a(xo, xi) = 7x0X1 

b = b{xo, xi) = SxqXi + + I j xo + Pxi, 

a,P,'y,5 E Zg, and 7 and 5 are selected such that $(a) t^t $(6) < 16. However notice that 
the coset representatives of the constructed cosets of RMg(l,m) generally correspond to cubic 
forms. 

It is noteworthy that a computer-based search failed to find kernels that generate more cosets 
of RMq(l,m) with PMEPR less than 4. Our search included all kernel functions in k variables 
for A; = 4 and q = 2, k = 3 and g = 4, and quadratic functions with k = 5 and q = 2. We 
leave the identification of more cosets of RMg(l,m) with low PMEPR (preferably close to 2) 
to further work. 

VIL Concluding Remarks 

In this paper we have shown how cosets of RMq(l, m) with low PMEPR can be obtained in 
a systematic way once suitable kernels are known. The PMEPR of these cosets is low because 
such cosets are comprised of sequences lying in so-called near-complementary pairs. It was 
demonstrated that those cosets of RMg(l,m) comprised entirely of Golay sequences simply 
arise as a special case in a more general theory. We presented some suitable kernels and, in this 
way, several previously unexplained phenomena arising in earlier works by Davis and Jedwab 
[4], by Paterson [23], and by Parker and Tellambura [19], [20] can now be understood in a 
general framework. Moreover we provided at least a partial answer to the question stated in 
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[23]: What other regularities appear in the PMEPRs of cosets as we move to higher alphabets, 
and how can they be explained in general? We have also established a connection between the 
OFDM and MC-CDMA coding problems, and suggested (generally nonbinary) coding solutions 
for MC-CDMA with low PAPR. 

We wish to point out some relations to the work by Parker and Tellambura [19], [20]. Indeed 
Theorem |7] in the present paper is similar to [20, Theorem 5] and [19, Theorem 6], which also 
exploit the Rudin-Shapiro construction. However these references do not identify the crucial 
connection between the PMEPR and the aperiodic auto-correlation, which allows us to extend 
the PMEPR bound to a slightly weaker bound in Theorem [TO] and Corollary [15] that holds 
for complete cosets of RM^(l,m). In contrast, in order to obtain cosets of RM2(l,m) with 
low PMEPR, the approach in [19] and [20] involves a computational search for the maximum 
PMEPR over a number of kernels. Indeed in [20, Table 6] it is implicitly shown that the cosets 
with coset representatives given in (|20l) have PMEPR at most 4 when q = 2. Notice that their 
semi-computationally obtained upper bound is exactly as predicted by Corollary \T5\ 

We presented some kernels of short length and we also sought good kernels of medium length. 
From the complexity point of view it is also realistic to perform an exhaustive search for good 
kernels of larger lengths, say 32 or 64, and defined over larger alphabets. We have not attempted 
such a search and leave it to further work. It would be desirable to find kernels producing cosets 
of RMq(l,m) with PMEPR close to 2, or even with PMEPR of exactly 2 (which implies the 
discovery of new complementary pairs). Moreover our theory would benefit from having an 
efficient way to construct good kernels. 
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